Genome-wide identification and expression analysis of ARF gene family in embryonic development of Korean pine (Pinus koraiensis).

BACKGROUND: The Auxin Responsive Factor (ARF) family plays a crucial role in mediating auxin signal transduction and is vital for plant growth and development. However, the function of ARF genes in Korean pine (Pinus koraiensis), a conifer species of significant economic value, remains unclear. RESULTS: This study utilized the whole genome of Korean pine to conduct bioinformatics analysis, resulting in the identification of 13 ARF genes. A phylogenetic analysis revealed that these 13 PkorARF genes can be classified into 4 subfamilies, indicating the presence of conserved structural characteristics within each subfamily. Protein interaction prediction indicated that Pkor01G00962.1 and Pkor07G00704.1 may have a significant role in regulating plant growth and development as core components of the PkorARFs family. Additionally, the analysis of RNA-seq and RT-qPCR expression patterns suggested that PkorARF genes play a crucial role in the development process of Korean pine. CONCLUSION: Pkor01G00962.1 and Pkor07G00704.1, which are core genes of the PkorARFs family, play a potentially crucial role in regulating the fertilization and developmental process of Korean pine. This study provides a valuable reference for investigating the molecular mechanism of embryonic development in Korean pine and establishes a foundation for cultivating high-quality Korean pine.

Comparative analysis of differential gene expression reveals novel insights into the heteroblastic foliage functional traits of Pinus massoniana seedlings.

Pinus massoniana needles, rich in medicinal polysaccharides and flavonoids, undergo heteroblastic foliage, transitioning from primary needles (PN) to secondary needles (SN) during growth, resulting in altered functional traits. Despite its significance, the molecular regulatory mechanisms governing these traits remain unclear. This study employs Iso-Seq and RNA-Seq analyses to explore differentially expressed genes (DEGs) associated with functional traits throughout the main growth season of heteroblastic foliage. Co-expression network analysis identified 34 hub genes and 17 key transcription factors (TFs) influencing light-harvesting antenna, photosystem I and II, crucial in photosynthesis regulation. Additionally, 14 genes involved in polysaccharide metabolism pathways, synthesizing sucrose, glucose, UDP sugars, and xylan, along with four genes in flavonoid biosynthesis pathways, regulating p-coumaroyl-CoA, quercetin, galangin, and myricetin production, exhibited differential expression between PN and SN. Further analysis unveils a highly interconnected network among these genes, forming a pivotal cascade of TFs and DEGs. Therefore, heteroblastic changes significantly impact needle functional traits, potentially affecting the pharmacological properties of PN and SN. Thus, these genomic insights into understanding the molecular-level differences of heteroblastic foliage, thereby establishing a foundation for advancements in the pharmaceutical industry related to needle-derived products.

Phytophthora pluvialis maintenance, spore production and detached needle assays.

Phytophthora pluvialis is an oomycete that primarily infects Pinus radiata and Pseudotsuga menziesii causing the destructive foliar disease red needle cast (RNC). Recent observations show that P. pluvialis can also infect western hemlock inducing resinous cankers. High-throughput and reproducible infection assays are integral to find key information on tree health and oomycete pathogenicity. In this protocol, we describe the propagation and spore induction of P. pluvialis, followed by detached needle assays for verification and quantification of virulence of P. pluvialis in P. radiata needles. These needle assays can be employed for high-throughput screening of tree needles with diverse genetic backgrounds. In downstream analysis, Quantitative PCR (qPCR) was utilized to assess relative gene expression, as exemplified by candidate RxLR effector protein PpR01. Additional techniques like RNA sequencing, metabolomics, and proteomics can be combined with needle assays and can offer comprehensive insights into P. pluvialis infection mechanisms.

Physiological responses to drought stress of three pine species and comparative transcriptome analysis of Pinus yunnanensis var. pygmaea.

Drought stress can significantly affect plant growth, development, and yield. Fewer comparative studies have been conducted between different species of pines, particularly involving Pinus yunnanensis var. pygmaea (P. pygmaea). In this study, the physiological indices, photosynthetic pigment and related antioxidant enzyme changes in needles from P. pygmaea, P. elliottii and P. massoniana under drought at 0, 7, 14, 21, 28 and 35 d, as well as 7 days after rehydration, were measured. The PacBio single-molecule real-time (SMRT) and Illumina RNA sequencing were used to uncover the gene expression differences in P. pygmaea under drought and rehydration conditions. The results showed that the total antioxidant capacity (TAOC) of P. pygmaea was significantly higher than P. massoniana and P. elliottii. TAOC showed a continuous increase trend across all species. Soluble sugar (SS), starch content and non-structural carbohydrate (NSC) of all three pines displayed a "W" pattern, declining initially, increasing, and then decreasing again. P. pygmaea exhibits stronger drought tolerance and greater recovery ability under prolonged drought conditions. Through the PacBio SMRT-seq, a total of 50,979 high-quality transcripts were generated, and 6,521 SSR and 5,561 long non-coding RNAs (LncRNAs) were identified. A total of 2310, 1849, 5271, 5947, 7710, and 6854 differentially expressed genes (DEGs) were identified compared to the control (Pp0D) in six pair-wise comparisons of treatment versus control. bHLH, NAC, ERF, MYB_related, C3H transcription factors (TFs) play an important role in drought tolerance of P. pygmaea. KEGG enrichment analysis and Gene set enrichment analysis (GSEA) analysis showed that P. pygmaea may respond to drought by enhancing metabolic processes such as ABA signaling pathway, alpha-linolenic acid. Weighted gene co-expression network analysis (WGCNA) revealed GST, CAT, LEC14B, SEC23 were associated with antioxidant enzyme activity and TAOC. This study provides a basis for further research on drought tolerance differences among coniferous species.

Transcriptome Study of Bursaphelenchus xylophilus Treated with Fomepizole Reveals a Serine/Threonine-Protein Phosphatase Gene that Is Substantially Linked with Vitality and Pathogenicity.

Bursaphelenchus xylophilus, the pine wood nematode (PWN), is the causal agent of pine wilt disease (PWD), which causes enormous economic loss annually. According to our previous research, fomepizole, as a selective inhibitor of PWN alcohol dehydrogenase (ADH), has the potential to be a preferable lead compound for developing novel nematicides. However, the underlying molecular mechanism is still unclear. The result of molecular docking showed that the stronger interactions between fomepizole and PWN ADH at the active site of ADH were attributed to hydrogen bonds. Low-dose fomepizole had a substantial negative impact on the egg hatchability, development, oviposition, and lifespan of PWN. Transcriptome analysis indicated that 2,124 upregulated genes and 490 downregulated genes in fomepizole-treated PWN were obtained. Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed genes indicated that fomepizole could be involved in controlling PWN vitality mainly by regulating key signaling pathways, such as the ribosome, hippo signaling pathway, and lysosome. Remarkably, the results of RNA interference indicated that the downregulated serine/threonine-protein phosphatase gene (stpp) could reduce the egg hatchability, development, oviposition, and lifespan of PWN, which was closely similar to the consequences of nematodes with low-dose fomepizole treatment. In addition, the silencing of stpp resulted in weakness of PWN pathogenicity, which indicated that stpp could be a potential drug target to control PWN.

Genetic parameter changes and age-age correlations in Pinus koraiensis growth over 40-year progeny testing.

BACKGROUND: Early selection in tree breeding could be achieved by addressing the longevity of tree improvement activities. Genetic parameter changes and age-age correlations are essential for determining the optimal timing of early selection. Practical tracking of genetic parameters of Pinus koraiensis, a major timber species with economic and ecological value, has become feasible as its progeny testing has entered the mid-term age in Korea. However, research on the age-age correlation of P. koraiensis as progeny trials approach rotation age is limited. This study aimed to investigate genetic parameter trends and age-age correlations in P. koraiensis progeny. P. koraiensis progeny were assessed at two sites using a linear mixed-effects model with two-dimensional spatial autoregressive structure. Height, diameter, and volume growth were measured in 11 assessments over 40 years. RESULTS: Genetic parameters, such as height and diameter, showed different patterns of change. The heritability ranged for the three growth traits in 0.083-0.710, 0.288-0.781, and 0.299-0.755 across the sites and age. Height heritability and its coefficient of variance decreased, whereas the diameter and volume estimates remained relatively constant. Correlations with Age 40 for phenotypic, genetic, and rank of breeding values ranged between 0.16 and 0.92, 0.594 and 0.988, and 0.412 and 0.965, respectively. These correlations generally increased as the age approached Age 40, with particularly high levels observed at Age 26 and Age 30. CONCLUSION: The observed genetic trends in P. koraiensis progeny testing offer valuable insights for early and precise selection. Notably, selecting superior genotypes at Ages 26-30 is supported by discernible genetic gains and robust correlations. Future research should integrate unbalanced data for selecting mother trees or families and conduct a comprehensive economic analysis of early selection to validate its practical benefits.

Response mechanism of carbon metabolism of Pinus massoniana to gradient high temperature and drought stress.

BACKGROUND: The carbon metabolism pathway is of paramount importance for the growth and development of plants, exerting a pivotal regulatory role in stress responses. The exacerbation of drought impacts on the plant carbon cycle due to global warming necessitates comprehensive investigation into the response mechanisms of Masson Pine (Pinus massoniana Lamb.), an exemplary pioneer drought-tolerant tree, thereby establishing a foundation for predicting future forest ecosystem responses to climate change. RESULTS: The seedlings of Masson Pine were utilized as experimental materials in this study, and the transcriptome, metabolome, and photosynthesis were assessed under varying temperatures and drought intensities. The findings demonstrated that the impact of high temperature and drought on the photosynthetic rate and transpiration rate of Masson Pine seedlings was more pronounced compared to individual stressors. The analysis of transcriptome data revealed that the carbon metabolic pathways of Masson Pine seedlings were significantly influenced by high temperature and drought co-stress, with a particular impact on genes involved in starch and sucrose metabolism. The metabolome analysis revealed that only trehalose and Galactose 1-phosphate were specifically associated with the starch and sucrose metabolic pathways. Furthermore, the trehalose metabolic heat map was constructed by integrating metabolome and transcriptome data, revealing a significant increase in trehalose levels across all three comparison groups. Additionally, the PmTPS1, PmTPS5, and PmTPPD genes were identified as key regulatory genes governing trehalose accumulation. CONCLUSIONS: The combined effects of high temperature and drought on photosynthetic rate, transpiration rate, transcriptome, and metabolome were more pronounced than those induced by either high temperature or drought alone. Starch and sucrose metabolism emerged as the pivotal carbon metabolic pathways in response to high temperature and drought stress in Masson pine. Trehalose along with PmTPS1, PmTPS5, and PmTPPD genes played crucial roles as metabolites and key regulators within the starch and sucrose metabolism.

Identification and Expression Patterns of WOX Transcription Factors under Abiotic Stresses in Pinus massoniana.

WUSCHEL-related homeobox (WOX) transcription factors (TFs) play a crucial role in regulating plant development and responding to various abiotic stresses. However, the members and functions of WOX proteins in Pinus massoniana remain unclear. In this study, a total of 11 WOX genes were identified, and bioinformatics methods were used for preliminary identification and analysis. The phylogenetic tree revealed that most PmWOXs were distributed in ancient and WUS clades, with only one member found in the intermediate clade. We selected four highly conserved WOX genes within plants for further expression analysis. These genes exhibited expressions across almost all tissues, while PmWOX2, PmWOX3, and PmWOX4 showed high expression levels in the callus, suggesting their potential involvement in specific functions during callus development. Expression patterns under different abiotic stresses indicated that PmWOXs could participate in resisting multiple stresses in P. massoniana. The identification and preliminary analysis of PmWOXs lay the foundation for further research on analyzing the resistance molecular mechanism of P. massoniana to abiotic stresses.

Genomics of the expanding pine pathogen Lecanosticta acicola reveals patterns of ongoing genetic admixture.

Lecanosticta acicola is the causal agent for brown spot needle blight that affects pine trees across the northern hemisphere. Based on marker genes and microsatellite data, two distinct lineages have been identified that were introduced into Europe on two separate occasions. Despite their overall distinct geographic distribution, they have been found to coexist in regions of northern Spain and France. Here, we present the first genome-wide study of Lecanosticta acicola, including assembly of the reference genome and a population genomics analysis of 70 natural isolates from northern Spain. We show that most of the isolates belong to the southern lineage but show signs of introgression with northern lineage isolates, indicating mating between the two lineages. We also identify phenotypic differences between the two lineages based on the activity profiles of 20 enzymes, with introgressed strains being more phenotypically similar to members of the southern lineage. In conclusion, we show undergoing genetic admixture between the two main lineages of L. acicola in a region of recent expansion. IMPORTANCE: Lecanosticta acicola is a fungal pathogen causing severe defoliation, growth reduction, and even death in more than 70 conifer species. Despite the increasing incidence of this species, little is known about its population dynamics. Two divergent lineages have been described that have now been found together in regions of France and Spain, but it is unknown how these mixed populations evolve. Here we present the first reference genome for this important plant pathogenic fungi and use it to study the population genomics of 70 isolates from an affected forest in the north of Spain. We find signs of introgression between the two main lineages, indicating that active mating is occurring in this region which could propitiate the appearance of novel traits in this species. We also study the phenotypic differences across this population based on enzymatic activities on 20 compounds.

Pine wilt disease: what do we know from proteomics?

Pine wilt disease (PWD) is a devastating forest disease caused by the pinewood nematode (PWN), Bursaphelenchus xylophilus, a migratory endoparasite that infects several coniferous species. During the last 20 years, advances have been made for understanding the molecular bases of PWN-host trees interactions. Major advances emerged from transcriptomic and genomic studies, which revealed some unique features related to PWN pathogenicity and constituted fundamental data that allowed the development of postgenomic studies. Here we review the proteomic approaches that were applied to study PWD and integrated the current knowledge on the molecular basis of the PWN pathogenicity. Proteomics has been useful for understanding cellular activities and protein functions involved in PWN-host trees interactions, shedding light into the mechanisms associated with PWN pathogenicity and being promising tools to better clarify host trees PWN resistance/susceptibility.

The impact of insect egg deposition on Pinus sylvestris transcriptomic and phytohormonal responses to larval herbivory.

Plants can improve their resistance to feeding damage by insects if they have perceived insect egg deposition prior to larval feeding. Molecular analyses of these egg-mediated defence mechanisms have until now focused on angiosperm species. It is unknown how the transcriptome of a gymnosperm species responds to insect eggs and subsequent larval feeding. Scots pine (Pinus sylvestris L.) is known to improve its defences against larvae of the herbivorous sawfly Diprion pini L. if it has previously received sawfly eggs. Here, we analysed the transcriptomic and phytohormonal responses of Scots pine needles to D. pini eggs (E-pine), larval feeding (F-pine) and to both eggs and larval feeding (EF-pine). Pine showed strong transcriptomic responses to sawfly eggs and-as expected-to larval feeding. Many egg-responsive genes were also differentially expressed in response to feeding damage, and these genes play an important role in biological processes related to cell wall modification, cell death and jasmonic acid signalling. EF-pine showed fewer transcriptomic changes than F-pine, whereas EF-treated angiosperm species studied so far showed more transcriptional changes to the initial phase of larval feeding than only feeding-damaged F-angiosperms. However, as with responses of EF-angiosperms, EF-pine showed higher salicylic acid concentrations than F-pine. Based on the considerable overlap of the transcriptomes of E- and F-pine, we suggest that the weaker transcriptomic response of EF-pine than F-pine to larval feeding damage is compensated by the strong, egg-induced response, which might result in maintained pine defences against larval feeding.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers.

MAIN CONCLUSION: The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Like mother like son: Transgenerational memory and cross-tolerance from drought to heat stress are identified in chloroplast proteome and seed provisioning in Pinus radiata.

How different stressors impact plant health and memory when they are imposed in different generations in wild ecosystems is still scarce. Here, we address how different environments shape heritable memory for the next generation in seeds and seedlings of Pinus radiata, a long-lived species with economic interest. The performance of the seedlings belonging to two wild clonal subpopulations (optimal fertirrigation vs. slightly stressful conditions) was tested under heat stress through physiological profiling and comparative time-series chloroplast proteomics. In addition, we explored the seeds conducting a physiological characterization and targeted transcriptomic profiling in both subpopulations. Our results showed differential responses between them, evidencing a cross-stress transgenerational memory. Seedlings belonging to the stressed subpopulation retained key proteins related to Photosystem II, chloroplast-to-nucleus signalling and osmoprotection which helped to overcome the applied heat stress. The seeds also showed a differential gene expression profile for targeted genes and microRNAs, as well as an increased content of starch and secondary metabolites, molecules which showed potential interest as biomarkers for early selection of primed plants. Thus, these finds not only delve into transgenerational cross-stress memory in trees, but also provide a new biotechnological tool for forest design.

Mitochondrial genome provides species-specific targets for the rapid detection of early invasive populations of Hylurgus ligniperda in China.

BACKGROUND: Hylurgus ligniperda, a major international forestry quarantine pest, was recently found to have invaded and posed a serious threat to the Pinus forests of the Jiaodong Peninsula in China. Continuous monitoring and vigilance of the early population is imperative, and rapid molecular detection technology is urgently needed. We focused on developing a single-gene-based species-specific PCR (SS-PCR) method. RESULTS: We sequenced and assembled the mitochondrial genome of H. ligniperda to identify suitable target genes. We identified three closely related species for detecting the specificity of SS-PCR through phylogenetic analysis based on 13 protein-coding genes (PCGs). Subsequently, we analyzed the evolution of 13 PCGs and selected four mitochondrial genes to represent slow-evolving gene (COI) and faster-evolving genes (e.g. ND2, ND4, and ND5), respectively. We developed four species-specific primers targeting COI, ND2, ND4, and ND5 to rapidly identify H. ligniperda. The results showed that the four species-specific primers exhibited excellent specificity and sensitivity in the PCR assays, with consistent performance across a broader range of species. This method demonstrates the ability to identify beetles promptly, even during their larval stage. The entire detection process can be completed within 2-3 h. CONCLUSIONS: This method is suitable for large-scale species detection in laboratory settings. Moreover, the selection of target genes in the SS-PCR method is not affected by the evolutionary rate. SS-PCR can be widely implemented at port and forestry workstations, significantly enhancing early management strategies and quarantine measures against H. ligniperda. This approach will help prevent the spread of the pest and effectively preserve the resources of Chinese pine forests.

Stochastic processes and changes in evolutionary rate are associated with diversification in a lineage of tropical hard pines (Pinus).

The study of the patterns of polymorphism and molecular evolution among closely related species is key to understanding the evolutionary forces involved in the diversification of lineages. This point is a big challenge in species with slow evolutionary rates, long life cycles, and ancient, shared polymorphisms such as conifers. Under the premise of divergence in a stepwise migration process, we expect clinal geographical patterns of purifying selection efficiency, and genetic structure related to latitude or longitude. If migration is accompanied by changes in the environment, we could further expect a role of positive selection in driving species divergence. Here, we infer patterns of polymorphism, efficiency of purifying selection, and molecular evolution using a dataset of 161 nuclear genes (∼71 Kb) in a lineage of hard pines from North America, the Caribbean, Mexico, and Central America presumed to have migrated from North America toward lower latitudes with tropical conditions. Under the premise of differences in selective pressures, we also look for possible signals of positive selection. To test our hypothesis, first we estimated different indices to infer patterns of polymorphism and efficiency of purifying selection (Ka, Ks, Ka/Ks, dN, dS, dN/dS, and dxy) and compared these metrics across five clades. Also, we investigated possible clinal patterns in these indices and morphological traits (needle length and cone length). Then we inferred genetic structure and environmental differences among species to test for possible signals of positive selection using phylogenetic methods in specific clades. We found differences among clades using Ka, Ks, and Ka/Ks with a relaxation of purifying selection, especially in the Elliotti and Patula clades. We also found environmental differences related to geographic distance, and among clades suggesting differences in selective pressures. The indices Ks, dxy, and needle length had relationships with geography but not ovulate cone length. Finally, we found that most analyzed genes are under purifying selection, but there was an exception of faster evolutionary rate in some pine species, suggesting the possible action of positive selection in divergence. Our study indicated that stochastic processes have played a key role in the diversification of the group, with a possible input of positive selection in pines from Mexico and Central America.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development.

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Effects of exogenous GA3 on stem secondary growth of Pinus massoniana seedlings.

Gibberellins (GAs) play a crucial role in regulating secondary growth in angiosperms, but their effects on the secondary growth of gymnosperms are rarely reported. In this study, we administered exogenous GA3 to two-year-old P. massoniana seedlings, and examined its effects on anatomical structure, physiological and biochemical changes, and gene expression in stems. The results showed that exogenous GA3 could enhance xylem development in P. massoniana by promoting cell division. The content of endogenous hormone (including auxins, brassinosteroids, and gibberellins) were changed and the genes related to phytohormone biosynthesis and signaling pathway, such as GID1, DELLA, TIR1, ARF, SAUR, CPD, BR6ox1, and CYCD3, were differentially expressed under GA3 treatment. Furthermore, GA3 and BR (brassinosteroid) might act synergistically in promoting secondary growth in P. massoniana. Additionally, lignin content was significantly increased after GA3 treatment accompanied by the express of lignin biosynthesis related genes. PmCAD (TRINITY_DN142116_c0_g1), a crucial gene involved in the lignin biosynthesis, was cloned and overexpressed in Nicotiana benthamiana, significantly promoting the xylem development and enhancing stem lignification. It was regarded as a key candidate gene for improving stem growth of P. massoniana. The findings of this study have demonstrated the impact of GA3 treatment on secondary growth of stems in P. massoniana, providing a foundation for understanding the molecular regulatory mechanism of stem secondary growth in Pinaceae seedlings and offering theoretical guidance for cultivating new germplasm with enhanced growth and yield.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

The nematicide emamectin benzoate increases ROS accumulation in Pinus massoniana and poison Monochamus alternatus.

Pine wilt disease (PWD) is caused by the pine wood nematode (PWN, Bursaphelenchus xylophilus) and transmitted by a vector insect, the Monochamus alternatus. The PWN has caused much extensive damage to pine-dominated forest ecosystems. Trunk injection of emamectin benzoate (EB) has been found to be the most useful protective measure against the PWN, due to its low effective dose and long residence time in the field. However, the interactions between EB and the host or the environment remain largely unknown, which limits the efficacy and stability of EB in practical field settings. In this study, we investigated the impact on PWN from EB injection for both adult and young host plants (Pinus massoniana) by taking a multi-omics (phenomics, transcriptomics, microbiome, and metabolomics) approach. We found that EB injection can significantly reduce the amount of PWN in both living adult and young pine trees. Additionally, EB was able to activate the genetic response of P. massoniana against PWN, promotes P. massoniana growth and development and resistance to Pine wilt disease, which requires the presence of PWN. Further, the presence of EB greatly increased the accumulation of reactive oxygen species (ROS) in the host plant in a PWN-dependent manner, possibly by affecting ROS-related microbes and metabolites. Moreover, we uncovered the function of EB limiting the consumption of P. massoniana by the JPS. Based on biochemical and gut microbial data, we found that EB can significantly reduces cellulase activity in JPS, whose transcription factors, sugar metabolism, and the phosphotransferase system are also affected. These results document the impact of EB on the entire PWD transmission chain through multi-omics regarding the dominant pine (P. massoniana) in China and provide a novel perspective for controlling PWD outbreaks in the field.

Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress.

Pinus massoniana is an important coniferous tree species for barren mountain afforestation with enormous ecological and economic significance. It has strong adaptability to the environment. TEOSINTE BRANCHED 1/CYCLOIDEA/PCF (TCP) transcription factors (TFs) play crucial roles in plant stress response, hormone signal transduction, and development processes. At present, TCP TFs have been widely studied in multiple plant species, but research in P. massoniana has not been carried out. In this study, 13 PmTCP TFs were identified from the transcriptomes of P. massoniana. The phylogenetic results revealed that these PmTCP members were divided into two categories: Class I and Class II. Each PmTCP TF contained a conserved TCP domain, and the conserved motif types and numbers were similar in the same subgroup. According to the transcriptional profiling analysis under drought stress conditions, it was found that seven PmTCP genes responded to drought treatment to varying degrees. The qRT-PCR results showed that the majority of PmTCP genes were significantly expressed in the needles and may play a role in the developmental stage. Meanwhile, the PmTCPs could respond to several stresses and hormone treatments at different levels, which may be important for stress resistance. In addition, PmTCP7 and PmTCP12 were nuclear localization proteins, and PmTCP7 was a transcriptional suppressor. These results will help to explore the regulatory factors related to the growth and development of P. massoniana, enhance its stress resistance, and lay the foundation for further exploration of the physiological effects on PmTCPs.